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Abstract 

We briefly comment on a paper by Rubano and Scudellaro [ astro-ph/0103335 
where they found general exact solutions for two classes of exponential po¬ 
tentials in a scalar field model for quintessence. In that paper the authors 
were led to some interesting conclusions after a proper choice of the integra¬ 
tion constants. By using dimensionless variables we show that the integra¬ 
tion constants can be found explicitly without additional assumptions. In 
consequence we revise some results and conclusions in that paper. We also 

reproduce observations for Type-la supernovae with good accuracy. 
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I. INTRODUCTION 


In a recent paper Rubano and Scudellaro found general exact solutions for two classes 
of exponential potentials in a scalar field model for quintessence. In that paper the au¬ 
thors studied a two-component perfect fluid (dust plus a scalar held) and considered an 
exponential potential and a combination of two exponential potentials for the quintessence 
held respectively. The quintessense action is given by S' = / d!^x + C^ + Cm}, 

where Cm is the Lagrangian for the matter degrees of freedom and the Lagrangian for the 
quintessense held is given by C^ = ~ ^(0)- Two classes of exponential potentials 

are of interest: 

Ui(0) = (1.1) 

and, 

1/2(0) = ( 1 . 2 ) 

where cx^ = and and R^ are generic constants. 

In Ref. [H] a hat (k=0) FRW universe was studied and the corresponding held equations 
that are derivable from an action principle are the following, 

2-= -/(/}. (1.4) 

Qj Qj o ^ 

and 

0 + 3-0-I-IT'(0) = 0, (1.5) 

a 

where D is the amount of matter (see Ref. [Q), the dot means derivative in respet to the 
cosmic time, the comma denotes derivative in respect to 0, and j = 1, 2. In that paper 
Rubano and Scudellaro were able to hnd general exact solutions to the above system of 
equations by introducing a pair of new variables: 
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3 J. ^ T / ^ N 

a =uv, (p = -inf —j, 

cr V 


in the first case (potential Vi, Eq. (1)), and 


= 


V? — v‘^ 


1 B{u + v)- 
, (p = - ln[- 


4 ’ ^ a y4(M — n) ’ 

in the second case (potential V 2 , Eq. (2)). The general solutions they found are: 

_2 jd2 —2 d 2 

U[t)=Uit + U2, v[t) = ——Uit H- —U2t +Vit + V2, 

0 2 

for the hrst class of potential Vi{(f)) and, 


( 1 . 6 ) 


(1.7) 


( 1 . 8 ) 


u{t) = + /de , v{t) = Vi sin {ujt + V 2 ) = hi sin (cut) + V 2 cos (cut), (1.9) 


for a combination of two exponentials V 2 ( 0 ), where uj = ABa^. 

In Ref. [0] the authors studied different situations by properly choosing the integration 
constants in Eqs. (1-8) and (1.9). The point of the present paper is to show that the 
integration constants in Eqs. (1.8) and (1.9) can be explicitly found without making addi¬ 
tional assumptions if we introduce dimensionless variables. In effect, let us introduce the 
dimensionless time variable r = where t is the cosmic time and Hq is the present 
value of the Hubble parameter, and the dimensionless scale factor a(r) = In these 

variables H{t) = a(r)/a(r) where, from now on, the dot means derivative in respect to the 
dimensionless time r. Then we have that, at present (r = 0), 


a(0) = 1, a(0) = 1 ^ H{0) = 1. 


( 1 . 10 ) 


Besides the changes t ^ ^ and a(t) a(0) a(r) in equations (l)-(5), one should 
replace ^ ^ - 3 ^ = D, where is the density of matter. 

After this rescaling, one can check that a^D = and, besides, Friedmann Eq. (1.3) 

evaluated at r = 0 implies for quintessence = 1 — VLrao ■ The net result is that equations 
(1.1)-(1.5) are unchanged but the constants and parameters, for instance uj = ABa^, are 
now dimensionless. While Ending the four integration constants in Eqs. (1.8) and (1.9) 
respectively, we can use the two conditions in Eq. (1.10) and two of the field equations 
(1.3)-(1.5) (for instance (1.3) and (1.4)) evaluated at r = 0. 
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II. AN EXPONENTIAL POTENTIAL 


For the potential of the class Vi (Eq. (1.1)), the solutions are 


u{t) = UiT + M2, v{t) = ^UiT^ + + ViT + V 2 , 

0 2 


The integration constants are the following 


(±) _^v/3 2{2-qo)-m 


=±— 


mo 


a‘^B^ 


( 2 . 1 ) 


( 2 . 2 ) 


,(±) 


.(±) 


V 2 --(±y> “i[±]- 

U 2 


3 + [±]^^2(l + go)-30 


mo (±) 

U 2 , 


(2.3) 


and, 


.,(±) _ 

[±] - 


2uf> 


(2.4) 


where go = ~(1 H{0)) is the present value of the deceleration parameter and, the (±) and 

[±], allow, in principle, for four different branches of the solution (2.1). 

The relation between the cosmological time f of M and our conformal time r is 


, 1 . M 2 ^ 

t = ■7ri'^ + —, 

Ho ui 


( 2 . 6 ) 


Substituting the above result in equation (1.8) and comparing with our solution gives the 
relation between our more general integration constants and those (denoted with a bar here) 
found in the particular case of [jl[: 


TT - - TT / - tt2 

Ui = HoUi,U 2 = 0,Vi = Ho{vi - -—),r;2 = f2 - 7 —= ujHq 


( 2 . 6 ) 


2 ui 6 mi 

Due to signs combination, indeed there are only two different solutions for the scale 


factor: 


a^(r) = 


{2>±RiYRl 4 {2>±Ri)Rl 3 2RI-RI + Q 
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H + 


-H + 3r + 1, 


(2.7) 


where 
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( 2 . 8 ) 


y/^\j2{l + go) ~ 3f2mo) — qo) — Silmo) 

Ri = -^-, R2 = -^-, 

For Ri to be real, in Eq. (2.8), Y^^2(T+'go)'^-'3i2^ should be real, then the following 
constraint holds. 


Qo ^ —1 + 


3. 

2^ 


mo • 


(2.9) 


In the limiting situation in which go = = 0? which implies that survive 

just one branch of the solution. 

In general, from Eqs. (2.8) and (2.7), one sees that the relevant magnitudes characterising 
the evolution of the universe depend only on two parameters Qq and flmo- These are not 
sensible to the particular value the dimensionless constant a; = a^B‘^ takes. 

Now we check the singular points of the solutions found. The condition a(r) = 0 means 
u{t) = 0 and/or v{t) = 0, i. e., we have four roots in r. The one derived from the condition 
u{t) = 0 is 


•in 777 j - 

3 + [±]f^2{l + qo)-3nmo 

If Tin were positive then we were faced with a very serious problem since r = 0 at present. 
However a careful analysis of (2.10) shows that, for the positive branch, is always negative, 
meanwhile for the negative branch, it could be positive only if the following condition holds: 
go > 5 + 3ffmg/2. But this condition is incompatible with present accelerated expansion. 
In case f (r) = 0 we found three roots, two of which are imaginary and one real. The real 
root is very complicated and we will not to write it here. A careful analysis shows that this 
root for Tin is negative as well, besides, by its absolute magnitude it is always greater than 
Tin in Eq. (2.10) (i. e., it is prior to Tin in Eq. (2.10)). Besides during the period between 


^This condition can be written, also, in the equivalent form go = Hqq, where the quintessence 

field stands for a dynamical cosmological constant. This relationship is often used in the bibliog¬ 
raphy on accelerated expansion 
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both roots the scale factor becomes imaginary so we drop times prior to Tin in (2.10). This 
is evident from Figs. 1 and 2, where the behavior of the scale factor is shown for both 
branches of the solution for the particular values hlmo = 0.3 and go = —0.4 and go = —0.45 
respectively. Similar arguments hold for the double exponential potential. The behavior of 
the scale factor of the universe shows that a(r) has zeroes for positive r in no case. 

Now we proceed to analyze how our solution reproduces experimental results. With this 
purpose, in Fig. 3 we plot the distance modulus S{z) vs redshift z, calculated by us and the 
one obtained with the usual model with a constant A term. The relative deviations are of 
about 0.5%. 


III. A MORE COMPLICATED POTENTIAL 


For the second potential V2{4>) (Eq. (1.2)) we write solutions in the following way: 

u{t) = ae‘^'^ + , n(r) = ui sin (cur) + ^2 cos (cur), (3.1) 


Here the situation is more complex. There are 8 branches of the solution. We studied 
one of these branches, leaving for further investigation a detailed study of the remaining 
ones. In the case studied here the integration constants a, f3, Ui, and V 2 in the solution (3.1) 
are found to be 


Vl = 


Y^4 + ^2^2(9 — 911^0 — 3(1 — 2go)) + 6^2 


4cu 


V2 = 


'3(1-2go) + 9-911, 


mo 


4cu 


a = JA + vi - (3, (3 = 


cu(4 + V 2 ) — 6 — UJV 1 V 2 


2ct^y 4 -f- V2 


(3.2) 


As before the situation simplihes if we consider the limiting case when go = —1 + |llmo- 


In this case the above integration constants look like 


3^'2 


ii-n 


mo 


2cu 


(3.3) 


We have three parameters in the model with this potential (A, B and cu). For simplicity, we 
have set cu = 1 in the following calculations. 
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Figure 4 shows the behaviour of the scale factor for flmo = 0.3 and go = —0.45, while 
hgure 5 plots 6{z) vs redshift z and compares with the usual experimental model with a 
constant A term |^. Relative deviations in this case are of about 0.25%. 

IV. CONCLUSIONS 

We showed that using the technique of adimensional variables and shifting the time 
origin, it was possible to hnd more general solutions than those found in Ref [|I|. 

In most cases we have used the simplifying condition go = —1 + which is justihed 

at least by the fact that a good agreement with experimental results is achieved. Indeed, we 
detected very small variations of 6{z) with relatively wide variations of ffmo; always within 
the aforementioned condition. This led us to explore another measurable quantities like the 
age of the universe. Equation (2.5) may be writen 

Hot = T + ^. (4.1) 

Ml 

Evaluating the above equation for the present (r = 0), 

Hoto = —, (4.2) 

Ml 

where to is the age of the universe. Using the equations (2.2) and (2.3) and the simplifying 
condition go = —1 + fflmo yields 

Hoto = -, (4.3) 

This is in agreement with the accepted standard cosmological paradigm {Hoio should be 
of order unity). However, better agreement is obtained if we use instead the more general 
requirement (2.9). As an illustration, consider the solution (2.7) for the single exponetial 
potential. From Eq. (2.10) it is seen that, for the negative branch of the solution (2.7) 
I < Hoto < 1 in all cases. For the positive branch Hoto < | and so, from experiment we see 
that the negative branch of the solution is more appropiated than the positive one. 


7 


As seen from the figures, the model for quintessence studied in Ref. [|l| yields an eternally 
accelerating universe with an event horizon that seems to be incompatible with superstring 
theory [^. One possible way to make this model compatible with observational evidence for 
a presently accelerating universe and with the absence of event horizons, is to add a negative 
constant term to the potentials Vj, equivalent to having a negative cosmological constant 
i- In a forthcoming paper we explore this possibillity. 

This short paper is dedicated to the memory of our colleague and friend Angelo Gino 
Agnese who introduced us to this subject. 

We acknowledge Andro Gonzales and Livan Rivero for help in computations and useful 
comments. We thank the MES of Cuba by hnancial support of this research. 
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FIGURES 

FIG. 1. The evolution of the scale factor is shown for the two branches of the solution for the 
single exponential potential for Qmo = 0-3 and qo = —0.4. 

FIG. 2. The evolution of the scale factor is shown for the two branches of the solution for the 
single exponential potential for Umg = 0.3 and go = —0.45. The two roots for which a(r) = 0 are 
visible. For the period between both roots the scale factor is imaginary, so this region is physically 
meaningless 

FIG. 3. The modulus distance 5{z) is ploted as a function of the redshift z. Dots represent the 
experimental values according to the standard metodology for Umo = 0.3. The solid line represents 
the theoretical curve obtained in the present model if we consider a single exponential potential 
for Qmo = 0.3 and go = —0.55 

FIG. 4. The evolution of the scale factor is shown for the two branches of the solution for the 
double exponential potential for Qmo = 0.3 and go = —0.45. 

FIG. 5. The modulus distance 5{z) is ploted as a function of the redshift z. Dots represent the 
experimental values according to the standard metodology. The solid line represents the theoretical 
curve obtained in the present model if we consider a single exponential potential. In both cases 
Umo = 0.3. We chose go = —0.55 
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